Disorder effects on resonant tunneling transport in GaAs/(Ga,Mn) As heterostructures 



(N 

o 

< 



Christian Ertler* 1 and Walter Potz 1 

institute of Theoretical Physics, Karl-Franzens University Graz, Universitatsplatz 5, 8010 Graz, Austria 

Recent experiments on resonant tunneling structures comprising (Ga,Mn) As quantum wells [Ohya 
et al., Nature Physics 7, 342 (2011)] have evoked a strong debate regarding their interpretation as 
resonant tunneling features and the near absences of ferromagnetic order observed in these struc- 
tures. Here, we present a related theoretical study of a GaAs/(Ga,Mn)As double barrier structure 
based on a Green's function approach, studying the self-consistent interplay between ferromagnetic 
order, structural defects (disorder), and the hole tunnel current under conditions similar to those in 
experiment. We show that disorder has a strong influence on the current-voltage characteristics in 
efficiently reducing or even washing out negative differential conductance, offering an explanation 
for the experimental results. We find that for the Be lead doping levels used in experiment the 
resulting spin density polarization in the quantum well is too small to produce a sizable exchange 
splitting. 

PACS numbers: 85.75.Mm, 73. 23. Ad, 73.63.-b, 72.25.Dc 
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I. INTRODUCTION 

Dilute magnetic semiconductors (DMS) are produced 
by doping of semiconductors with transition metal el- 
ements, which provide local magnetic moments aris- 
ing from open electronic d or / shells. 0, Bulk 
Gai-^Mn^As may be regarded as the prototype: Mn re- 
siding on the Ga site (Mnc a ) donates a hole, associated 
with valence band p-orbitals, and provides a local mag- 
netic moment associated with partly filled Mn d-orbitals. 
MnQa is a moderately deep acceptor with the energy lev- 
els lying about 100 meV above the valence band edge.[|| 
By increasing the Mn density the acceptor levels become 
more and more broadened, developing into an impurity 
band which allows hole propagation and, for sufficiently 
high doping level, is believed to merge with the valence 
band. 0] At the same time, Mn more and more takes un- 
wanted lattice positions, such as the antisite and inter- 
stitial position in the fee lattice, or may even form Mn 
clusters, all leading to strong electron-hole compensation 
which eventually destroys ferromagnetic ordering. There 
is some debate as to the order in which these events occur 
as the Mn concentration is increased. Probably due to 
the presence of unintentional defects in (Ga,Mn)As sam- 
ples, depending on growth conditions, experimental evi- 
dence has led to somewhat conflicting conclusions about 
the precise position of the Fermi level in ferromagnetic 
bulk (Ga,Mn)As.@] Some experiments can be interpreted 
by placing it into the top of a GaAs-like valence band 
edge which is broadened by disorder. [l| Others suggest 
the existence of an isolated impurity band in the ferro- 
magnetic state. 0-0] 

Recently a systematic series of experiments in form 
of non-equilibrium tunneling spectroscopy on double- 
barrier resonant tunneling structures with a (Ga,Mn)As 
quantum well 0-Q were conducted to provide a deeper 
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insight into this question. The group reported a near ab- 
sence of ferromagnetic order in the well under bias and 
obtained weak signatures of resonant tunneling, observ- 
able only in the second derivative of the current-voltage 
(IV) characteristic. Their conclusion that the Fermi en- 
ergy lies in the impurity band has evoked strong debates 
and an alternative explanation has been given, which 
proposes that the resonant-tunneling signature is caused 
merely by the confined states in a potential pouch formed 
at the contact/barrier heterointerface. jlOj In this expla- 
nation the observed dependence of the peak positions 
on the quantum well width is completely attributed to 
the increased series resistance which, however, seems to 
be insufficient to account for all well-width-dependent 
trends in the experimental results, as discussed in detail 
in a reply by Tanaka et al. which again emphasizes the 
existence of quantized levels in the (Ga,Mn)As quantum 
wells, [m 

Indeed, quantization effects can be expected in 
(Ga,Mn) As for a layer thickness of about < 3 mn since in 
a recent scanning tunneling microscopy experiment the 
radius of the Mn acceptor wave function has been de- 
termined to be about 2 nm H and one can expect that 
near the band edge Bloch-like and delocalized eigenstates 
will coexist in the picture of merging impurity and va- 
lence bands. [l2l| Tunneling spectroscopic experiments of 
(Ga,Mn)As quantum well structures have indicated such 
effects. [7| However, the signatures in the current -voltage 
characteristics appear to be rather weak and no regions of 
negative differential resistance due to resonances associ- 
ated with (Ga,Mn)As well layers have been observed as of 
yet, with the notable exception of an asymmetric magne- 
toresistance resonant-tunneling structure. 13] This sug- 
gests that a significant concentration of unwanted defects 
and/or disorder may be present, depending on growth 
conditions, as it is known to be the case in thin layers of 
amorphous Si, in which similarly weak signatures have 
been found. [Til [lj| The density of imperfections due to 
the presence of Mn interstitial or antisite defects can be 
as high as 20% of the nominal Mn doping, which makes 
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(Ga,Mn)As a heavily compensated system. [l6|, Oj} Even 
lower structural quality for (Ga,Mn)As must be expected 
at heterointerfaces since the need of low-temperature epi- 
taxy for growing the (Ga,Mn)As layers is harmful to 
forming clean interfaces with other materials. Moreover, 
the interstitial defects may be trapped near the inter- 
faces in post-growth annealing procedures which have 
been found successful for bulk (Ga,Mn)As. This suggests 
that transport through thin layers of (Ga,Mn)As is in- 
fluenced by disorder and defects more severely than in 
annealed bulk structures. 

The growth of heterostructures, on the other hand, 
provides the appealing opportunity to drive (Ga,Mn)As 
layers into a genuine non-equilibrium situation by means 
of an external bias which modifies their local hole den- 
sity, pos sibly leading to bias-dependent ferromagnetic 
behavior. [l8l [l9| However, drawing conclusions from the 
physics of a thin (Ga,Mn)As layer regarding the Fermi 
energy position in the bulk is an intricate problem, since 
a Fermi energy in a (Ga,Mn)As quantum well under bias 
conditions is not well defined. 

In a recent series of studies we have investigated the 
ferromagnetic bias anomaly in (Ga,Mn)As-based het- 
erostructures and reached the conclusion that, for suf- 
ficiently high hole densities in the thin (Ga,Mn)As quan- 
tum wells, ferromagnetic ordering becomes bias depen- 
dent leading to variable spin-polarized currents. |l8l lilij 
Here, we study the low doping regime (relative to the 
Mn concentration) and use a refined model for the valence 
band states which accounts for both heavy and light -hole 
states. This allows us a direct comparison to recent ex- 
periments and, as will be shown, enhances the effect of 
disorder on suppressing a resonant-tunneling signature 
in the IV characteristic. Based on a four band Kohn- 
Luttinger Hamiltonian the transport properties are inves- 
tigated within a self-consistent non-equilibrium Green's 
function method which accounts for space charge effects 
and a hole-density-dcpendent exchange splitting. We 
show that disorder reduces or even completely washes 
out regions of negative differential conductance in the IV 
curve. We find that, for the Be lead doping levels as 
used in experiment, the resulting spin density polariza- 
tion in the quantum well is low and thus leads to almost 
vanishing ferromagnetic order. Our theoretical model is 
presented in Sect. |TT] and the results and relevance to 
experiment are discussed in Sect. IIIII Summary and con- 
clusions are given in Sect. flV] 



II. PHYSICAL MODEL 

Here we describe our transport model for heterostruc- 
tures composed of layers of GaAs, GaAlAs, and 
(Ga,Mn)As grown along the z-axis. In this study the 
band structure of the top of the valence bands is modeled 
by the Kohn-Luttinger Hamiltonian [201 ] , which allows us 
to take into account the mixing of heavy hole (HH) and 
light hole (LH) bands, which is of crucial importance for 



getting a realistic transmission function for holes tun- 
neling through a double-barrier structure as shown in 
Ref. [2l|. Ordering the four spin-3/2 basis vectors at 
the T-point as u a = \^,m a ) with m a = {l)4> — \i~ fL 
the wave-vector-dependent Kohn-Luttinger Hamiltonian 
reads 
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The matrix elements can be expressed in terms of the 
dimensionless Luttinger parameters 71,72 and 73: 
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where m is the free electron mass. In order to consider- 
ably simplify the numerical demands for the calculation 
of macroscopic quantities, such as the current density, 
which require the summation over the in-plane momen- 
tum, we apply the axial approximation in which the con- 
stant energy surface in the k-space becomes cylindrically 
symmetric but for which HH-LH band mixing is still in- 
cluded. Within the axial approximation the transmis- 
sion function only depends on the absolute value of the 
in-plane momentum k 2 p = k x + ky. Space-dependent (in 
z-direction) potentials are taken into account within the 
envelope function approximation, which effectively leads 
to replacing k z — » by —id z . By approximating the intro- 
duced spatial derivatives on a finite grid of spacing a one 
ends up with an effective nearest-neighbor tight-binding 
Hamiltonian of tridiagonal form 

H =Yj e *l' c l,<? c l,<r' + c l+i,* c l,o' + h - c -> ( 3 ) 
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with c, denoting the creation operator for site / and 
orbital a. The on-site and hopping matrices, respectively, 
take the form 
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Here, the matrix elements are given by 
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with 7 = (72 + 7 3 )/2. This effective tight-binding model 
has the advantage that space-dependent potentials, ex- 
change splittings in magnetic layers, and structural im- 
perfections can be readily included in the orbital onsite 
energies of the model, i.e., the diagonal elements of the 
onsite matrix using 

= E** + Ul - - I A, + £ ra nd (7) 

with Ui denoting the intrinsic hole band profile due to 
the band offset between different materials, <f> is the the 
electrostatic potential, e is the elementary charge, A/ de- 
notes the local exchange splitting in the magnetic mate- 
rials with a = ±1, and e ran d labels a random shift due to 
disorder, as will be detailed below. 

With the ferromagnetic order being mediated by the 
itinerant carriers the exchange splitting of the hole bands 
self-consistently depends on the local spin density of the 
holes. It can be derived within an effective mean-field 
model taking into account two correlated mean magnetic 
fields stemming from the ions' d-electrons spin polariza- 
tion (Sg) and the hole spin density (s z ) = (n^— n\ )/2.p3- 
l24j ] The exchange splitting of the hole bands is then given 

by 

A(z) = -J p dn imp (z)(S z )(z) , (8) 

with z being the longitudinal (growth) direction of the 
structure, J p d > is the exchange coupling between the 
p-like holes and the d-like impurity electrons, and n- lmv (z) 
is the impurity density profile of magnetically active ions. 
The effective impurity spin polarization (S z ) is induced 
by the magnetic field caused by the mean hole spin po- 
larization, yielding 

(S.) = SBs (^^) , 0) 

where, ks is the Boltzmann constant, T is the lattice 
temperature, and B$ is the Brillouin function of order S, 
here with S = 5/2 for the Mn impurity spin. Combining 
the last two expressions gives the desired result 

A(.) = j^(z)SB 8 {^^M^ . (10) 



Since the hole spin density (s z ) is changed by the in- and 
out-tunneling holes, the magnetic and transport proper- 
ties of the system are coupled self-consistently. 

To obtain realistic potential drops between the two 
leads space-charge effects have to be taken into account. 
In the Hartree approximation the electric potential is de- 
termined by the Poisson equation, 

dd^ 
az az 

where e and N a , respectively, denote the dielectric con- 
stant and the acceptor density. The local hole density at 
site \l) can be obtained from the non-equilibrium "lesser" 
Green's function G < : 

"(0 = ^E / ^G<{E;laM) , (12) 

k||,<7 

with A and ky , respectively, being the in-plane cross sec- 
tional area of the structure and the in-plane momentum. 
The lesser Green's function is determined by the equation 
of motion 

G< = G R T,<G A (13) 

where G R and G A = [G fl ] + denotes the retarded and 
advanced Green's function, respectively. The scattering 
function E< = S^+S^ describes the particle inflow from 
the left (I) and right (r) reservoirs [25| with 

E,< =/o(£- W ,r)(££ P -E*.) , (14) 

where fo(x) = [1 + cxp(x/fcsT)]~ 1 is the Fermi dis- 
tribution function and jii and /j, r , respectively, denote 
the quasi-Fermi energies in the contacts. The retarded 
and advanced self-energy terms = + E^ and 
E' 4 = [E fi ] + , respectively, couple the simulated system 
region to the left and right contacts. The surface Green's 
function of the leads is needed to obtain the contact self- 
energy T, R and is calculated by using the algorithm of 
Lopez-Sancho et al.[26| The retarded Green's function of 
the system, finally, is given by 

G R = [E + iri-Hs-^y 1 . (15) 

which we calculate by consecutively adding one layer of 
the system after another which, in our case, solely re- 
quires the inversion of a 4x4 matrix for each additional 
layer. 

The transport equations Eqs. (fT3|) and (| 1 5(1 couple via 
the spin-resolved hole density to the exchange splitting 
of the hole bands Eq. (|10p . and the Poisson equation 
Eq. (fTTj) . For a given applied voltage this system of equa- 
tions is solved in a self-consistent loop until convergence 
of the electrostatic potential and the exchange field is 
reached. A small external magnetic field is applied ini- 
tially to aid spontaneous symmetry breaking. For the 
next bias iteration, the self-consistent solution from the 
previous bias value is used for an initial guess. Having 



obtained the self-consistent potential profile the trans- 
mission probability T(E, kii) from the left to the right 
reservoir is calculated by 



log(LDOS) at V= 



T(£,k||) =Ti[TiG R r r G A ] 
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with Tr[-] denoting the trace operation and the lead cou- 
pling functions being defined by ^(i, r ) — i\^fi r ) ~ r )]- 
The steady-state current density is obtained by an in- 
tegration over all incoming k|| states and the total energy 
E: 

3 = YhJ2 I dETiE^lME^ri-ME,^)} (17) 



with the applied bias V = (/i; — /i r )/e being defined as 
the difference in quasi-Fermi levels of the contacts. 

We model disorder effects in the (Ga,Mn)As layers by 
performing a configurational average over structures with 
randomly chosen diagonal elements of the onsite matrix 
resulting in an ensemble of tight-binding Hamiltonians. 
The diagonal onsite matrix elements are sampled accord- 
ing to a Gaussian distribution around e ran d = in Eq. J?]) 
for increasing standard deviations (<7dis = 20, 40, 60, and 
80 meV). For each specific structure (and Hamiltonian) 
the current-voltage (IV) curve is solved self-consistently 
for an upsweep of the applied voltage. This effective one- 
dimensional modeling of disorder must be viewed as a 
limited estimate since it corresponds to a cross-sectional 
average of transport though uncorrelated effective lin- 
ear chains. As such, any disorder correlations parallel to 
the interface are neglected. Such correlations in disor- 
der will play a role in the establishing of ferromagnetic 
order in real structures relative to the idealized homoge- 
neous mean-field model adopted here, since both ferro- 
magnetic order and disorder effects are highly dependent 
upon spatial dimensionality. [13, HH Also, this type of av- 
eraging cannot model the effects of Mn clusters. Addi- 
tional types of disorder from Mn clustering, Mn intersti- 
tials, etc. however may be present in real structures, but 
their type and concentration may differ from sample to 
sample. 



III. 



RESULTS 



We investigate a double barrier structure consisting of 
Be-doped GaAs leads and a (Ga,Mn)As quantum well, 
similar to the experimental setups presented in RefsIH, |9|. 
The structural similarity of the layer materials allows 
one to use the same valence band model for the whole 
structure, which considerably simplifies the theoretical 
description. In recent resonant tunneling spectroscopic 
experiments of the group of Tanaka 0], however, a Schot- 
tky barrier contact with Au at one side and a GaAs:Be 
contact with an AlAs barrier on the other side of the 
(Ga,Mn)As quantum well was used. If the transport is 
primarily determined by the resonant levels in the well, 




FIG. 1: (Color online) Logarithm of the local density of states 
(LDOS) at zero bias for the perfect structure (no disorder) 
using an inverted energy scale. The self-consistent potential 
profile is indicated by the (red) solid line, whereas the Fermi 
energy position of the leads is plotted as (red) dashed line. 
The position of the impurity band in the well is indicated 
schematically by the (white) bold solid line below the valence 
band edge. 



our results obtained for a completely zincblende struc- 
ture will be relevant also for the Schottky-barrier system 
investigated in experiment. 

For the simulations we assume following parameters, 
comparable to the experimental values of Ref. H|: barrier 
thickness db&r = 4 layers (sa 2 nm), quantum well width 
d w — 10 layers (« 5 nm), barrier height Vbar = 300 meV, 
relative permittivity of GaAs e r = f2.9, exchange cou- 
pling constant J p d = 0.06 eV nm 3 , 5% Mn doping, and 
the lead temperature T = 4.2 K. The Fermi level in the 
GaAs leads is chosen to /!/ = 10 meV, which corresponds 
to a Be-doping of about f 18 cm~ 3 as used in experiment. 
The dimensionless Luttinger parameters are set to stan- 
dard values for GaAs: 71 = 6.85,72 = 2.f, and 73 = 2.9. 

While (Ga,Mn)As inevitably is a disordered system, 
our starting point is an idealized system for which the 
valence band edge is identical to that of GaAs. Sub- 
sequently, valence band disorder is added in increasing 
steps allowing for a systematic qualitative account of its 
consequences on the IV-curve. For convenience we use an 
inverted hole energy scale but retain positively charged 
holes. In order to simulate the presence of MnQ a impu- 
rity band levels, partially populated by holes, we use a 
positive (repulsive) charge background of 10 18 cm -3 in 
the (Ga,Mn) As layer leading to an upward shift of the po- 
tential profile in the well region, as shown in Fig. [TJ The 
local density of states of the ideal double-barrier struc- 
ture (in absence of disorder) at zero bias and V ~ 0.1 V, 
corresponding to the dominant current peak of the IV 
curve (see Fig. [5]), is shown in Fig. [T] and Fig. [3J respec- 
tively. The resonant levels in the well are clearly visible 
and the solid line indicates the self-consistent potential 
profile of the structure. At zero bias the valence band 
edge of (Ga,Mn)As is lifted relative to that of the contact 
GaAs layers by about 30 meV thus partially aligning the 
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FIG. 2: (Color online) Logarithm of the local density of states 
(LDOS) at the bias V = 0.1 V for the perfect structure (no 
disorder) using an inverted energy scale. The self-consistent 
potential profile is indicated by the (red) solid line. The Fermi 
energy position of the left and right lead is plotted as (red) 
dashed and dotted line, respectively. Holes of the impurity 
band (indicated schematically by the bold (white) solid line 
below the valence band edge) can tunnel out to the collector 
reservoir. 
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FIG. 3: (Color online) The transmission function at zero bias 
for the perfect structure (solid line), for moderate disorder of 
<7dis = 20 meV (dashed line), and stronger disorder <7di s = 
60 meV (dotted line) by sampling over 3 x 10 4 configurations. 



MriQa levels (indicated schematically by the bold (white) 
solid line) with the chemical potential which, in turn, lies 
about 10 meV above the valence band edge of the GaAs 
leads. Therefore, when a bias greater than about 10 mV 
is applied the Mnca levels can no longer be filled reso- 
nantly from the emitter side and, beyond about 30 mV, 
no longer from either emitter or the collector. The latter 
situation is shown schematically in Fig. [2j This loss of 
holes from the MnQ a levels in the well region and the 
insufficient resupply of holes from the GaAs emitter into 
the resonant valence band levels lead to a breakdown of 
ferromagnetic order in the (Ga,Mn)As well under small 
bias, i.e. a zero bias anomaly. 

The transmission function and local density of states 
in the center of the Gai-^Mn^As quantum well region 



.10' 



3- 

CD 

c 

en 
O 
a 



10 



10" 




a.. =0 

— dis 




a.. = 20 meV 

- dis 




a.. = 60 meV 
dis 





0.05 0.1 0.15 
Energy (eV) 



0.2 



FIG. 4: (Color online) Local density of states (LDOS) in 
the middle of the quantum well at zero bias for the perfect 
structure (solid line) , for moderate disorder of <7dis = 20 meV 
(dashed line), and stronger disorder <Tdi s = 60 meV (dotted 
line) by sampling over 3 x 10 4 configurations. 
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FIG. 5: (Color online) IV-characteristic of a magnetic double 
barrier structure for increasing degree of disorder measured 
by the standard deviation <7di s of the Gaussian distribution of 
the randomly chosen diagonal onsite matrix elements £ ran d • 



is plotted for increasing degree of disorder in Fig. [3] and 
Fig. 21 respectively. Disorder leads to defect level in the 
band gap and a significant spectral broadening of the 
resonant levels associated with the lowest valence band 
resonances. Note also the strong increase of transmission 
probability in the low-bias regime from disorder, opened 
by resonances for tunneling. 

As a key difference compared to the bulk (Ga,Mn)As 
situation, the hole density inside a (Ga,Mn)As quantum 
well is established by forming a steady state situation 
with the external leads. Since the Be doping level in the 
GaAs contacts is significantly lower than the Mn con- 
centration this leads to a strongly reduced hole concen- 
tration (over bulk) in the (Ga,Mn)As wells under bias. 
Even under favorable bias conditions, in which resonant 
levels in the quantum well can be populated from the ex- 
ternal leads, the quantum well hole density remains on 
the order of the hole density in the GaAs leads (in our 
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case « 10 18 cm -3 ), which is at least two orders of magni- 
tude smaller than in bulk (Ga,Mn)As. Therefore, in our 
simulations we find practically vanishing ferromagnetism 
(exchange splitting) for all bias values. In previous stud- 
ies we have shown that for higher hole doping of the 
leads a bias-dependent ferromagnetic state with a maxi- 
mum exchange splitting of the order of tens of meV can 
be expected, being detectable by a significant spin po- 
larization of the collector current. [HI, In this earlier 
study we focused on the transport through the first HH- 
subband using a simpler effective mass band model. The 
picture that the quantum well hole density is determined 
by the coupling to the leads surely applies to the case 
when the impurity band merges with the valence band, 
leading at most, to a broadening and shift of the valence 
band edge. In the presence of an isolated impurity band 
loosely bound holes may exist (at least at low voltages) in 
the confined impurity bands lying energetically below the 
first valence 2D-subband (on the inverted energy scale). 

The main result of the paper is given in Fig. [5j which 
shows the IV-characteristics for increasing degree of dis- 
order. Using a parallelized code typically 240 configura- 
tions are used for each characteristic, which needs about 
4 days of computation on a 12-node Opteron server for 
a single IV-curve. For increasing disorder the resonances 
in the IV-curve become more and more broadened and 
start to overlap until they are almost washed out. Here 
a model which takes into account the HH-LH band mix- 
ing in the well is essential in order to see this effect. 
If ignored, the LH resonances would dominate the cur- 
rent magnitude compared to the HH peaks and regions of 
pronounced negative differential resistance would persist 
even for unphysically high degrees of disorder. For the 
low doping regime of the leads (relative to the Mn concen- 
tration in the quantum well) , as considered here, we find 
an almost vanishing ferromagnetism (A < 1CV 2 meV) in 
the well leading to a vanishing current spin polarization. 
This results also suggest that, for a given (small) degree 



of disorder, enhancing the confinement, e.g. by using 
a thinner quantum well and/or higher barriers, should 
lead to smaller overlap between the subbands, thus en- 
hancing resonant-tunneling features the IV-curve. Al- 
ternative external pressure may be used to enhance the 
splitting between the lowest HH and LH subbands. 



IV. CONCLUSIONS 

In summary, we have shown that two factors can be 
relevant to understand the experimental findings of weak 
resonant tunneling features and an absence of ferromag- 
netic order arising from thin (Ga,Mn)As quantum-well 
layers: 

(i) The presence of considerable disorder in thin layers 
of (Ga,Mn)As which conspires with HH-LH band mixing 
in the quantum well to efficiently weaken any signatures 
of resonant tunneling in the IV characteristics. 

(ii) The observed (almost) vanishing of ferromagnetic 
order in the (Ga,Mn)As quantum well can be understood 
by finding orders of magnitude lower hole densities in the 
well as compared to the bulk (Ga,Mn)As case of equal 
Mn concentration. 

We thus find the original interpretation of observing 
quantum size effects in a (Ga,Mn)As quantum well by 
the group of Tanaka 0-Q to be plausible and consistent 
with our numerical results. 
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